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We develop a new, dynamical field theory of isotropic randomly branched polymers, and we use 
this model in conjunction with the renormalization group (RG) to study several prominent problems 
in the physics of these polymers. Our model provides an alternative vantage point to understand the 
swollen phase via dimensional reduction. We reveal a hidden Becchi-Rouet-Stora (BRS) symmetry 
of the model that describes the collapse (#-)transition to compact polymer-conformations, and 
calculate the critical exponents to 2-loop order. It turns out that the long-standing 1-loop results 
for these exponents are not entirely correct. A runaway of the RG flow indicates that the so-called 
^'-transition could be a fluctuation induced first order transition. 



PACS numbers: 64.60.ae, 05.40.-a, 64.60.Ht, 64.60. Kw 

A single linear (non-branched) polymer in solution un- 
dergoes a second order phase transition from a swollen 
to a collapsed state when the solvent temperature sinks 
below the so-called 0-point. In the swollen phase, the 
polymer can be thought of as a self-avoiding walk, and 
its radius of gyration or Flory radius scales with monomer 
number N as R N ~ N" SAW [y SAW > 1/2). In the col- 
lapsed phase, the polymer assumes a compact globule- 
like conformation, and Rn ~ N x / d where d is the dimen- 
sionality of space. The understanding of this collapse 
transition as a critical phenomenon has advanced con- 
siderably over the years [l[ . 

In comparison, much less is known about the col- 
lapse transition of randomly branched polymers (RBPs) . 
There exists a number of numerical studies [Hj S El Hi B 0] 
that, taken together, indicate that the phase diagram is 
fairly complex including a line of collapse transitions that 
has qualitatively distinct parts. One part, called the 9- 
line, corresponds to continuous transitions with univer- 
sal critical exponents of swollen RBP configurations with 
mainly tree-like character to compact coil-like configura- 
tions. The other part of the transition line, called the 9'- 
line, corresponds to the collapse of foam- or sponge-like 
RBPs to vesicle- like compact structures. In 2 dimensions, 
one finds nonunivcrsal exponents if this transition is con- 
sidered as continuous Q ■ The two different parts of the 
collapse-transition line are separated by a multicritical 
point which belongs to the isotropic percolation univer- 
sality class. One of the open questions is the existence 
of a possible further transition line between the differ- 
ent configurations of collapsed RBPs. As far as theory is 
concerned, it is the swollen phase that is best understood 
mainly because the statistics of swollen RBPs can be for- 
mulated in terms of an asymmetric Potts model [1, [j| [T3| 
although Flory theory and real space renormalization 
[l2T | have also been successfully applied. The former ap- 
proach was used in particular to solve the field theoretic 
problem via a mapping of the relevant part of the asym- 
metric Potts model to the Yang-Lee edge problem using 



dimensional reduction [13j. In contrast, the collapse of 
RBPs has been much less studied, and the current un- 
derstanding mainly rests on the seminal field theoretic 
work of Lubensky and Isaacson (LI) Q and Harris and 
Lubensky However, it turns out that these papers, as 
far as they consider the collapse (0-)transition, contain a 
fundamental error in the renormalization procedure, and 
as a consequence the long-standing 1-loop results for the 
collapse transition are strictly speaking not correct. In 
addition, it is not clear to date whether the 9' transition 
is a second order transition or not. Therefore, we feel 
that the important RBP problem deserves a fresh look. 

In this paper we develop a new, dynamical field the- 
ory [l5| for RBPs based on a model for dynamical per- 
colation with a tricritical instability [T(| in the non- 
percolating phase whose very large clusters (lattice ani- 
mals), at critical values of the control parameters, have 
the same statistics as collapsing RBPs. We discuss the 
relation of our model to the asymmetric Potts model and 
carefully analyze its symmetries. In the swollen phase, 
the model has a high super-symmetry including trans- 
lation and rotation invariance in super-space and leads 
to the well known Parisi-Sourlas dimensional reduction 
fl3j | . At the collapse transition, super-rotation symmetry 
is broken, and we only have translation invariance in su- 
perspace, i.e., Becchi-Rouet-Stora (BRS) symmetry fl8| . 
We perform a 2-loop renormalization group (RG) calcu- 
lation, that corrects and extends the long standing LI 
results for the collapse transition. Furthermore, we show 
that the ^'-transition is characterized by a runaway of 
the RG flow which suggests that this transition is a fluc- 
tuation induced first order transition contrary to what 
has been assumed in recent numerical studies [2J. 

Our field theory (for background on field theory meth- 
ods, we refer to [TH, [2(|) is based on a generalization 
of the general epidemic process. For a related approach 
to the somewhat simpler problem of directed randomly 
branched polymers, see [2l| . The primary fields of our 
theory are the field of agents n(r, t) and the field of the in- 
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active debris m(r, t) = X J_ dt' n(r, t') which ultimately 
forms the polymer cluster. The minimal non-Markoffian 
Langevin equations describing the process are given by 

f 1 
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C(r, t)((r>, t>) = [X-'gnir, t)5(t - t') - fn(r, t)n{v\ t') 

xS{r-r'). (lb) 

The parameter r tunes the "distance" to the percolation 
threshold. Below this threshold r is positive. The term 
proportional to c describes the influence of the debris 
on diffusion. For the ordinary percolation problem, this 
term is irrelevant. As long as g 1 > 0, the second order 
term f'm 2 is irrelevant near the transition point and the 
process models ordinary percolation. We permit both 
signs of g' so that our model allows for a tricritical in- 
stability. Consequently we need the second order term 
/' > for stabilization purposes, i.e., to limit the density 
to finite values. The process is assumed to be locally ab- 
sorbing, and thus all terms in the noise-correlation func- 
tion contain at least one power of n. The first part of the 
noise correlation takes into account that the debris arises 
from spontaneous decay of agents, and thus g > 0. The 
term proportional to / > simulates the anticorrclating 
behavior of the noise in regions where debris has already 
been produced. 

Now, we refine these Langevin equations into a field 
theoretic model for RBPs. This procedure involves a 
number of nontrivial steps that we will briefly sketch in 
the following and that will be presented in detail else- 
where [lil ]. As the first step, we represent the Langevin 
equations as a stochastic response functional 



J = I d d x{x I dt 
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in the Ito-sense [15|, |22|, [2J, [2J . This functional has the 
benefit that it allows us to systematically calculate av- 
erages (• • ■ ) of all sorts of observables via functional in- 
tegration with weight exp[— Cf\. For studying polymers, 
we focus on a single cluster of a given size N which we 
assume to emanate from a small source of strength q at 
the origin r = at time t = 0. Then, the key quantity is 
the probability distribution for finding a cluster of mass 
N given by [T^ 



V(N) ~ q(5(N - M)n(0, 0)) , 



(3) 



where M = J d d rm(r, oo). V(N)/N is expected to be 
proportional to the partition sum for interacting lattice 
animals Q up to an non-universal exponential factor 
~ Pq if N becomes large. In actual calculations, the 
delta function appearing in averages like in Eq. {3| is 
hard to handle. This problem can be simplified by av- 
eraging over Laplace-transformed observables, which are 



function of a variable conjugate to N, say z, and ap- 
plying inverse Laplace transformation in the end. The 
switch to Laplace-transformed observables can be done 
in a pragmatic way by augmenting the original J with 
a term zM. and then working with the new functional 
J z = J+zM.. Because we are interested here only in the 
static properties of the final cluster after the epidemic has 
become extinct, we can greatly simplify the theory by fo- 
cusing on the frequency zero part of J z , that is taking the 
quasistatic limit [HI, OS [2{| m(r,t) — > moo(r) = itp{v), 
n(r,t) — » —iip(r). Taking this limit, one has to be care- 
ful to account for the causal ordering of fields that results 
from the Ito calculus. In diagrammatic perturbation the- 
ory, this means that one has to rule out diagrams with 
closed propagator loops. An elegant way to achieve this 
is to use so-called ghost fields whose sole purpose is to 
generate additional diagrams that cancel any diagrams 
with non-causal loops. Such a procedure does not change 
the physical content of the theory but simplifies calcula- 
tions and makes it easier to find higher symmetries. The 
required cancellations can be achieved [lj| by using D 
commuting (bosonic) fields Xi subject to the constraint 
12iLi Xi — so that they form the irreducible represen- 
tation (-D, 1) of the permutation-group Sd, and taking 
the limit D — * — 1 at the end of the calculation. Further- 
more, we eliminate redundant parameters by rescaling, 
mixing, and shifting the fields. After all, we obtain the 
quasistatic Hamiltonian 

+ i(roX (2) + (V X ) (2) ) + ! f^(2w+X (2) ) 
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where we use the shorthand notation = 12iLiXii 
and where h is a shifted version of the Laplace variable 
z. The r's and the g's are combinations of the original 
parameters, cf. Eqs. {1} and ([2]). In particular, tq and t\ 
are linearly related to r and g', respectively, so that, in 
mean-field theory, the collapse transition corresponds to 
vanishing To, t±, and h and swollen RBPs correspond to 
vanishing ro and h, and positive and finite n. 

What is the connection between our Hamiltonian (H|) 
and other, established models for RBPs, Percolation and 
the Yang-Lee problem? To address this question, we 
rescale the fields so that g[ = go, (which is possible, of 
course, only if both arc non-vanishing, in particular at 
RG fixed points), and we define a new order parameter 
field with (D + 2) components, s\ = i(p, S2 = itp, and 
for i > 3: Sj = Xi-2 — [si + S2)/D. This field satis- 
fies the Potts constraint Y^i=i s i = 0' an d the resulting 
Hamiltonian with Sd+i permutation-symmetry is that of 
the asymmetric {D + 2)-state Potts model which lies at 
the heart of the known formulations of the RBP prob- 
lem [t| [n|. For <?i = .92 = 0, the model reduces 
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to the symmetric (D + 2)-state Potts model with Sr>+2- 
symmctry and thus produces the field theory of percola- 
tion in the limit D — > — 1. For go + 2g\ = go + Ag 2 = 0, 
the Hamiltonian decomposes in a sum of D + l uncoupled 
Hamiltonians each describing the Yang-Lee edge prob- 
lem. 

To reveal the connection of our work to the results by 
Parisi and Sourlas [l3| for swollen RBPs and to shed light 
on the collapse transition from a symmetry perspective, 
it is interesting to discuss the super-symmetries of our 
model. If g[ is zero (or irrelevant like for finite n > 
0) non-causal loops are isolated, and can therefore be 
eliminated with a pair of Fermionic ghost fields ip and ip 
p^ . Using anticommuting super-coordinates 9, 9 with 
integration rules Jd91 = j d9 1_=0, / d6 9 = J d9 9 =1 
and defining a super-field <3>(r, 9,9) = ip(r) + 9ip(r) + 
-0(r)# + 99(p(r), we can recast our model Hamiltonian as 



+ i 



{i$(r -V 2 -Tid 9 -d e )$ + i2$ 
* 3 + Y^idgdg^) - ^{dgde®) 2 ] } . (5) 



This Hamiltonian shows BRS-symmetry [18j,l20|, i.e., 7i ss 
is invariant under a super-translation 9 —> 9 + e, 9 — > 
9+e. Moreover, if the control parameter t\ is positive and 
finite, i.e., if we consider the problem of swollen RBPs, t\ 
can be reset by a scale transformation to 2. The super- 
coordinates become massive, and the derivatives combine 
to a super-Laplace V 2 +Tid§de — > V 2 + 2d§de =: □. The 
coupling constants, g\ and g 2 become irrelevant, and can 
be neglected. Then the Hamiltonian takes the super- 
Yang-Lee form and attains, besides the super-translation 
invariance, super-rotation invariance. Now dimensional 
reduction [l3[ can be used to reduce the problem to the 
usual Yang-Lee problem in two lesser dimensions which 
culminates into to well known results for swollen RBPs. 

Now we come to the heart of our RG analysis, where we 
focus on the case that the control parameters tq and T\ 
take critical values (zero in mean-field theory) where the 
correlation length diverges, and correlations between dif- 
ferent polymers vanish. The actual objects of our pertur- 
bation theory are the naively UV-divergent vertex func- 
tions Tj^ , which consist of irreducible diagrams with k 
and k amputated legs of <p and ip, respectively, as func- 
tions of the wavevector q. We calculate these functions in 
dimensional rcgularization and e-expansion about d = 6 
dimensions (e = 6 — d) to 2-loop order and then remove 
their UV divergences in minimal subtraction using the 
scheme 

if, <P, X) -> X) = Z 1/2 (<p + Kip, ip, x) , (6a) 

Tj -> fi = Z~ 1 Z t3 T j , (6b) 

h^h = Z- 1 ' 2 (h + IfJ.-^TiAijTj) , (6c) 

g a ^ g a = Z- 3 / 2 (u a + B a )^ 2 , (6d) 

where [i is an inverse length scale that is used to make the 
coupling constants dimensionlcss, where (tj) = (tq,ti) 



and (g a ) = (go, g[, g±, g?). Note that the renormalization 
scheme introduces a counter term proportional K that 
has no counterpart in the Hamiltonian (j4|). This term 
can be viewed as a remnant of the term proportional to c 
in the original response functional ([2]) which we removed 
in our journey towards 7i because c is in the sense of the 
RG redundant. As a counter term this term is indispens- 
able, however, because the quadratically divergent ver- 
tex function L 2 ,o(q) =r2,o(0) + q 2 L 2 ' (0) + . . . contains 
a UV-divergent T'2 (0). This fact was overlooked by LI 
Q in their calculation, and their long-standing 1-loop re- 
sults are incorrect although, fortunately, the numeric de- 
viations from the correct 1-loop results are rather small. 

As it stands, the Hamiltonian (j4|) has a remaining 
rescaling invariance that makes one of the coupling con- 
stants redundant. Before we can analyze the RG flow, we 
need to remove this redundancy. To this end, we switch 
to rescaling invariant fields <p —* ff^ 1 ^, <p —* gofi, con- 
trol parameters r — ► r , T\ — * go 2 T\, h — ► (2g )^ 1 h, 
and effective dimensionlcss coupling constants u = woit'i, 
v = UqUi, w = UqU 2 - The fixed points of the RG flow 
are determined by the zeros of the Wilson functions for 
the three effective couplings, (3 U = [id^u\o (|o indicates 
that unrenormalizcd quantities are kept fixed while tak- 
ing derivatives) and so on. Our calculation produces 

Pu = (-£ + 7u/2 + 10v)u , (7a) 

(3 V = (-£ + 25u/6 + 2lv/2)v - 5w/6 , (7b) 
(i w = (-2e + 21u/2 + 2bv)w 

- (hu 2 + 2Quv/2+llv 2 )v , (7c) 

where we refrain from showing the 2-loop parts of our 
results due to space constraints. The picture of the RG 
flow that arises from these equations is the following: 
The BRS-planc u = is an invariant surface of the flow 
equations (|7all7cl) to all orders and divides the (u, v,w)- 
space in two parts: the percolation-part with u > 
and the Yang-Lee-part with u < which is non-physical 
for the branched polymer problem. The percolation line 
v = w = is an invariant line for both signs of u. For 
u > the flow goes to the percolation fixed point whereas 
for u < the flow tends to infinity. The Yang-Lce-line 
with a = b = 0, where a = u + 2v and b = it ; +4to, is also 
an invariant line for both signs of u. For u < the flow 
goes to the Yang-Lee fixed point whereas for u > the 
flow runs away to infinity. Altogether we have six fixed 
points which are compiled in Table [T] Besides the trivial 
Gaussian fixed point we find in the BRS-plane the stable 
collapse fixed point (named Collapse), and an instable 
fixed point (Inst2). This point lies on a separatrix in the 
BRS-plane and is attracting on it. The flow of the part 
which contains the collapse fixed point is of course at- 
tracting to Collapse. The other part shows runaway flow. 
Turning to the percolation-part of the (it, v, u>)-space, we 
find one instable point (Percolation) on the percolation 
line v = w = 0. Because Percolation has two stable di- 
rections it defines a plane that divides the space in two 
parts. The flow in one of it goes to Collapse whereas the 
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TABLE I: RG fixed points to leading order. 



flow in the other part is again running away. The sta- 
bility plane of Percolation for u > is a continuation of 
the separatrix found above on the BRS-planc for u = 0. 
In the Yang-Lee-part of the (u, v, u>)-space, we also find 
a plane which is the continuation of the BRS-separatrix 
now into the region with u < 0. This plane is separated 
in two parts by the Yang-Lce-line. One part is attracting 
to an instable fixed point (Instl), the other part shows 
runaway flow. Both planes divide the (u, v, w)-space in 
a wedge-shaped part attracting to Collapse, and a part 
where the flow goes to infinity. The edge of the wedge is 
the separatrix in the BRS-planc. Note that the flow di- 
agram has the following perhaps unexpected implication 
for the ^'-transition. The region behind the percolation 
plane where u runs away to ever more positive values 
indicate that this transition might be discontinuous and 
not, as previously assumed, a second order transition. 
Finally, we compile our main results for the collapse 



transition. Our RG analysis leads to three indepen- 
dent critical exponents. For the probability distribution 
V(N), we find the asymptotic form 



V(N,y)/N ~N- e ^f v (yN<t>), 



(8) 



where /iq is a non-universal constant, f-p is a scaling func- 
tion, and the effective control parameter for the "dis- 
tance" from the transition is given by the scaling vari- 
able y which is a linear combination of tq and t\. For 
the radius of gyration, we obtain 

R{N,y)~N v f R (yN*). (9) 

To second order in e-cxpansion, the critical exponents of 
the collapse transition read 



= - - 0.4925 (e/6) - 0.5778 (e/6) 2 , 



0.0225 (e/6) - 0.3580 (e/6) 2 , 



v =- + 0.1915 (e/6) + 0.0841 (e/6)' 



(10a) 
(10b) 
(10c) 



Note that these results compare well within the expecta- 
tions for large e with recent simulations in d = 2 [2j . 

In summary, we have presented a new renormalized 
field theory for RBPs. Though almost a classic physics 
problem, RBPs are still a lively subject of current re- 
search with important open questions some of which our 
work can help to settle. 
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